Bethe ansatz equations for Bariev's correlated electron chain with boundaries 
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The Bethe ansatz equations are presented for Bariev's correlated electron chain with boundaries. 
■ This is achieved by using the coordinate space Bethe ansatz method. 
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This work concerns with Bariev's correlated electron model jjj with boundary terms. The quantum integrability 
of this boundary model has been established by one of the present authors || using Sklyanin's formalism for the 
boundary quantum inverse scattering method ||. The aim of this report is to diagonalize this model by means of 
the coordinate space Bethe ansatz technique, which has been used to solve the Hubbard open chain jl]-[6) and the 
. supersymmetric U model with boundary terms ||]. We will derive the Bethe ansatz equations. Our results may 
be used to investigate the finite size corrections to the low-lying energies in the system. 

Our starting point is the Hamiltonian of Bariev's correlated electron model with boundaries: 
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^bulk _|_ ^boundary 

L-1 

^buik = ^[( c t T c j+1T + ct +1T c iT ) exp[T)nj +n ) + (ctjCj-m +c] +n Cj l )exp(T]n^)}, 

^boundary = Q+ exp ( 3?? ) [ 2 s i n h T 1 n N] n N[ + CX P (-7 ? ) (n N1 + n Nl )] 

H [2sinh?7ni T n u + exp(-ry)(ni T + n u )]. (1) 
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Here cj Q and c JQ are, respectively, the creation and annihilation operators of electrons with spin a(=f or J.) at site j 
and the number density operator. 

In order to derive the Bethe ansatz equations for the model, we assume that the eigenfunction of Hamiltonian (|l|) 
■ takes the form 

m 

O: l*)= £ *<n,-,vAxw-,XN)ti in ■■■ct N(TN \0), 

00 

; ^ai,-,aN(xi,---,x N ) = ^epA aQU ...^ QN {k P Qi, - ■ • ,k P Q N )e-x^(i^k Pj Xj), (2) 

ON '■ p i= l 

where the summation is taken over all permutations and negations of ki , • • • , kjy , and Q is the permutation of the N 
particles such that 1 < xqi < ■ ■ ■ < xqn < L. The symbol ep is a sign factor ±1 and changes its sign under each 
'mutation'. Substituting the wavefunction into the the eigenvalue equation H\^f) = E\^}, one may get 

£^ A. . . ,a j ,ai , • ■ ■ (■ ■ ■ j kj , ki , • • •) — Sij (ki , kj )A. . . ai a ^ _ . . . (• • • , ki, kj , • • • ) , 



A aii ...(-kj, ■■■) = s (k j ;p la JA a ^...(k J , ■ ■ ■), 



A.. )(Ti (- ■ ■ , -kj) = s R {k j ;p L<Ji )A...^(- ■ ■ , kj), (3) 
where Sij is the two-particle scattering matrix, and s L and s R are the boundary scattering matrices: 
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Sij (ki, kj 



q sinh(i 2 ) e' 2 sinh -q g 

sinh(i fcl 2 fc2 +77) sinh(i fcl 2 ^ 2 +7?) 

q e 2 sinln; smh(i 2 ) q 

sinh(i fcl 2 fc2 +7?) sinh(i kl 2 ^ 2 +7?) 
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In the above equations, 



P\a=P\=C_e 2n , pLa = PL — c + e 2v . (5) 

For later use, let us introduce the notation 



l—pe lk 
1 — pe l 



Then, the boundary scattering matrices s L and s R may be expressed as 

s L (fc i ;pi (Ti ) = s(-kj;p 1(Ti ), 

s R (k f ,p L(7i ) = s(k r ,p L(Ti )e ik ^ L+1 \ (7) 

As is seen from the above equations, the two-particle S matrix is nothing but the R-matrix of the asymmetric six- vertex 
model, which satisfies the quantum Yang-Baxter equation (QYBE), 

Sij(ki,kj)Sii(ki,ki)Sji(kj,k{) = Sji(kj,ki)Sn(ki,ki)Sij(ki,kj) (8) 

whereas the boundary scattering matrices s L and s R are propotional to unit matrices, which satisfy the reflection 
equations, 

Sji( — kj,—ki)s (kj) Plcr^Sij ( fej, kj)s (ki\P\<ri) 

= s {ki ; picn )Sji ( kj 7 ki)s (kj)pi C r i )Sij(ki ) kj) ) 
Sji (—kj,—ki)s (kj) PLuj ) Sij (ki, kj ) s [ki] PLai ) 

= s R (k i -p Lai )S ji (k J , -k i )s R (k j ;p LcT% );p crz )S ji (k j ,k l ). (9) 

Indeed, substituting (0) into ([}]), one obtains the consistency contition for the matrix s(kj)p cri ), 

Sij (ki, kj)s(ki] Pat) Sji (kj , —ki)s(kj] p ai ) 

= s(kj',p a ^)Sij(ki, kj)s(ki',p rTi )Sji( kj, ki). (10) 

This is just the reflection equation discussed in Mezincescu and Nepomechie's work || (see also Q), which admits 
a trivial solution s k 1. Then, the diagonalization of Hamiltonian ([j]) is reduced to solving the following matrix 
eigenvalue equation 

Tjt = t, j = l,---,N, (11) 
where t denotes an eigenvector on the space of the spin variables and Tj takes the form 

Tj = 57(%) S i (-%;^)i?7(fc j )flt(fc J 0^(fci;P^)5/(fc j ) (12) 

with 

Sj'ikj) = S j ^ N (k :j ,k N ) ■ ■ ■ Sjj + i(kj,k j+ i), 
Sj (^j) — Sj,j—x (kj , kj—i) ■ ■ ■ Sj : i(kj, k\), 
Rj(kj) = Sij(kx, -kj) ■ ■ ■ Sj-i t j(kj-i,-kj), 

R t( k j) = Sj+i, J (k ]+ i,-k ] ) ■ ■ ■ S N .j(k N , -kj). (13) 

This problem may be solved using the algebraic Bethe Ansatz method. Indeed, it can be shown that [Tj, Tj] =0 for 
any i, j = 1, • • • , N, which implies that Tj, j — 1, 2, • • • , TV are simultaneously diagonalizable. We will not present 
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the proof here, but mention that Tj is related to Sklyanin's transfer matrix t(A) for an inhomogeneous asymmetric 
six-vertex model with open boundary: 



e -ik 3 2(L + l) 

r(X)=tr (MT (X)To 1 (-X)), (14) 

where 

£( fc 5P) = \ ^ lk ' 
1 — pe lK 

T (A) = Sbjv(A - k N ) ■ ■ ■ 5 i(A - fci), 

sinh(iA + t?) / e-" \ . . 

M= sinh( i A + 2 ?7 ) { e"J- (15) 

Therefore, the commutativity of the Sklyanin's transfer matrix r(A) for different values of A guarantees the commu- 
tativity of Tj'a. Diagonalizing r(A) using the algebraic (or analytical) Bethe ansatz method and setting A = kj, we 
derive the following Bethe ansatz equations 



H sinh(i^^ + 2) sinh(i^4^ + 2) 



s=i sinh(i J — -) sinh(i 3 2 "° — -) 

* sin h(i^^ + |) S inh(z^±^ + j) _ TT sinh(z^^ + V ) sinh(i^±^ + rj) 

sinh(i^^ - 2) sinh ( z I^i -|) ~ J-l sinh^^p _ ,,) S inh(i^±^ -77) ' (16) 

With the solutions {fcj} of @, the energy eigenvalue E 1 of the model is given as 

AT 

£ = 2^cos%. (17) 
j'=i 

In summary, we have presented the Bethe ansatz equations for the ID Bariev's correlated electron chain with 
boundaries in the context of the coordinate space Bethe ansatz technique. Our results would be useful in analysing 
the structure of the ground state and some low lying excitations of the model in the thermodynamic limit. An 
interesting problem is to investigate spectrum of boundary states in the model, as done in ]lO[ for the open Hubbard 
chain. 
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